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We analyzed effects of elasticity on the dynamics of fluids in porous media by studying a flow 

■ ' " ' of a Maxwell fluid in a tube, which oscillates longitudinally and is subject to oscillatory pressure 

M^ \ gradient. The present study investigates novelties brought about into the classic Biot's theory 

^--s . of propagation of elastic waves in a fluid-saturated porous solid by inclusion of non-Newtonian 

fvj ' effects that are important, for example, for hydrocarbons. Using the time Fourier transform and 

, , transforming the problem into the frequency domain, we calculated: (A) the dynamic permeability 

'~{ • and (B) the function F[n) that measures the deviation from Poiseuille flow friction as a function of 

^~~i ' frequency parameter k. This provides a more complete theory of flow of Maxwell fluid through the 

y—^ , longitudinally oscillating cylindrical tube with the oscillating pressure gradient, which has important 

f^ ' practical applications. This study has clearly shown transition from dissipative to elastic regime in 

which sharp enhancements (resonances) of the flow are found. 
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bJO; I. INTRODUCTION 

. rH ! ^ quantitative theory of propagation of elastic waves in a fluid-saturated porous solid was formulated in the classic 

^ ■ papers by Biot [1]. One of the major findings of Biot's work was that there was a breakdown in Poisseuille flow above 
a certain characteristic frequency 

specific to the fluid-saturated porous material. Biot theoretically studied this phenomenon by considering the flow 

of a viscous fluid in a tube with longitudinally oscillating walls under an oscillatory pressure gradient. Apart from the 

fundamental interest, the investigation of the dynamics of fluid in porous media under oscillatory pressure gradient 

and oscillating pore walls is of prime importance for the recently emerged technology of acoustic stimulation of oil 

reservoirs [2]. For example, it is known that the natural pressure in an oil reservoir generally yields no more than 

approximately 10 percent oil recovery. The residual oil is difficult to produce due to its naturally low mobility, and the 

f^ ■ enhanced oil recovery operations are used to increase production. It has been experimentally proven that there is a 

t^^ . substantial increase in the net fluid flow through porous space if the latter is treated with elastic waves. However, there 

^^ ' is a fundamental lack of understanding of the physical mechanisms of fluid mobilization in saturated rock through the 

^ , effect of elastic waves; the theory of such mobilization virtually does not exist. Biot's theory can be used to describe 

:^ ■ the interaction of fluid-saturated 

c/3 \ solid with the sound for a classic Newtonian fluid; however, oil and other hydrocarbons exhibit signiflcant non- 

^Z^ ' Newtonian behavior [3]. The aim of this paper is therefore to incorporate non- Newtonian effects into the classical 
^; study of Biot [1]. 

Recently, del Rio, Lopez de Haro, Whitaker [4] presented a study of enhancement in the dynamic response of a 

viscoelastic (Maxwell) fluid flowing in a stationary (non-oscillating) tube under the effect of an oscillatory pressure 

gradient. We combine this theory with the effect of the acoustic oscillations of the walls of a tube introduced by Biot 

. !^ ■ [l]i providing a complete description of the interaction of Maxwell fluid, filling the pores, with acoustic waves. 

K> Finally, in order to emphasize that the consept of dynamic permeability is an adequate way of desciprtion of the 

; I phenomenon, we note that this concept has been widely used before [5]. 

C^ ' In the next section we formulate our model, whereas the last section will conclude with the discussion of the results. 

II. THE MODEL 

In this section we present our model of Maxwell fluid flowing in a cylindrical tube whose walls are oscillating 
longitudinally and the fluid is subject to an oscillatory pressure gradient. We give analitical solutions of the problem 
in the frequency domain. 

The governing equations of the problem consist of the continuity equation for the incompressible fluid, 

V-t/^O, (1) 



and the linearized nioinentuin equation, 



dv 



-Vp - Vf. 



(2) 



Here, v, p, p denote velocity, pressure and mass density of the fluid, whereas f represents the viscous stress tensor. 
We describe the viscoelatic effects of the fluid using Maxwell's model, which assumes that 



ot 



(3) 



where "q is the viscosity coefficient and tm is the relaxation time. 

Now, let w be a velocity of the wall of the tube which oscillates in time as e~*'^*, where w is the angular frequency. 
The flow of fluid in a cylindrical tube with longitudinally oscillating walls can be descibed by a single component of 
the velocity, namely, its z-componcnt v^ {z axis is along the cetcrlinc of the tube) . We use the cylindrical coordinate 
system (r, (/), z) in the treatment of the problem. We introduce the relative velocty Ui as Ui = v^—u. Thus, assuming 
that all physical quantities vary in time as 

e~*"*, we arrive at the following master equation for Ui 



V^Ui 



Ul = (1 - ILOtm). 

V V 



(4) 



Here, we have introduced the following notations: 



X 






which is a sum of the applied pressure gradient and force exerted on the fluid from the oscillating wall of the tube 
and v^ which Va v = r]/ p. Note that a no-slip boundary condition at the wall is assumed. 
The solution of Eq. (4) can be found to be [1] 

t/i = -- + CJo(/3r), 

lUJ 



where Jq is the Bessel function and j3 = \/{uP-tm + ioj)/v. 

Applying a no-slip boundary condition Ui{a) = at the wall of the tube, where a is its radius, we finally get 



Ui{r) 



X 
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MM 



MPa) 
Defining the cross-section averaged velocity as 
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Ui{r)rdr, 
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we obtain 



Ul 



Xa'^{l — ILUtjn) 1 



1 



2Ji(/3a) 
{pa)Jo{(3a) 



K{uj)X. 



(6) 



Here K{ijj) is the dynamic permeability [4] that describes the frequency dependent response of the tube to the applied 
total force on the fluid. Simple comparison reveals that Eq.(6) resembles closely Eq.(6) from Ref. [4], with the only 
difference being that we have X in place of their dp/dz. Note that in the case of stationary tube walls {p{du/dt) -^ 0) 
Eq.(6) exactly coincides with Eq.(6) from Ref. [4]. Simple calculation (applying L'Hospital rule for the 0/0 uncertanty) 
shows that linit^^o Kioj) — a? /{%v). Thus, following Ref.[4], we will introduce the dimensionless dynamic permeability 
as K*{uj) — K{llj)/K{0), which will be used later (see Fig. 7). Note, that we were able to easily reproduce Fig. 1 from 
Ref. [4], confirming the existence of sharp resonances of K*{lu) in the elastic regime (see below for the definition of 
this regime). 

Following the work of Biot [1] we calculate the stress at the wall r. 
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1 — iutr. 



dUijr) 
dr 
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iu}{l — itotm) Jo{(3a) 



(7) 



Note, that when tm —^ this expression obviously coincides with the corresponding Newtonian form. 

The total friction force is 27raT. Following Biot we calculate the ratio of total friction force to the average velocity, 
i.e. 



27rar 2TTf]{/3a)[Ji{/3a)/Ja{Pa)] 
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2Ji(/3a) 



Simple analysis reveals that 



(1 - iuitm) 



lim —^^ ~ OTTT], 
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(8) 
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which corresponds to the limiting case of Poiseuille flow. Following Biot [1], we also introduce a function F{k) with 
K being n ~ ay/ujjv in the following manner 



2'KaT 



= 87rr/F(K), 



thus. 



F{n) 



1 Ky^i + K-^ /a J\(K^i + K?- ja)! J^in^i + n^ ja) 



(1 — iti^ / a) 



2Ji{n^i + K^/a) 



Kyji + K^ jaJi^^K^i + K^ ja) 



(9) 



as 



Note, that F{n) measures the deviation from Poiseuille flow friction as a function of frequency parameter k, 
introduced by Biot [1]. 

In Eq. (9) , OL denotes the Deborah number [4] , which is defined as the ratio of characteristic time of viscous effects 
it, = (P' jv to the relaxation time t^, i.e. a = ty/tm = o?' /{vtm)- 

As noted in Ref. [4], the value of the parameter a determines in which regime the system resides. Beyond a certain 
critical value [ac = 11.64), the system is dissipative, and viscous effects dominate. On the other hand, for small a 
{a < ac) the system exhibits viscoelastic behavoir which we call the elastic regime. 

Note, that the Newtonian flow regime can be easily recovered from Eq.(9) by putting a ^ oo. We plot this limiting 
case in Fig. 1. As it can be seen from the figure, both Re[F{K)] and /to[F(k)] exactly coincide with the Newtonian 
limiting case studied in Blot's work (see Fig. 4 in Ref.[l]). This graph demonstrates a breakdown in Poiseuille fiow 
as frequency increases (recall that k oc ^/uj) ■ In all our calculations we have used polynomial expansions of Jq and Ji 
with absolute error not exceeding 10^^ %. Thus, our calculation results are accurate to this order. 

A finite-but-large a regime is shown in the next two figures. Fig. 2 corresponds to the case when a = 10*. We 
see in Fig. 2 that the real and imaginary parts of F(k) start to deviate from the Newtonian fluid behavior at large 
frequencies. In Fig. 3 solutions correspond to the case when a = 100.0; thus we see how viscoelastic effects become 
pronounced already at low frequencies. 

Fig. 4 presents the behavior of Re[F(K)] and Im[F{K)] when a — 10.0. As it can be seen from the graph, sharp 
resonances appear on the curves. This feature can be explained by the fact that in this case a is less than ac, which 
means that the system switched from the dissipative (viscous) regime to the elastic one. 

Fig. 5 shows the solutions when a = 1.0. In this case we see more irregular behavior of Re[F{K)] and Im[F{K)] 
with a number of irregular spikes. 

The extreme non-Newtonian (elastic) regime is studied in Fig. 6 where we plot the solutions for the case when 
a = 10"'^. In this case notable change is that there are fewer but more pronounced spikes. Re[F{K)] is close to zero 
for most of the frequencies, and only at certain frequencies we see sharp resonances. In this regime the system acts a 
window for these frequencies. 

Another noteworthy observation is that, on one hand, as long as a > ac = 11.64 (Figs. 2 and 3), Re[F{K)] is always 
greater than its initial value, i.e. Re[F{K)] > 1, and for large frequencies it reaches a ceratin asymptotic value. On 
the other hand, when a < ac (Figs. 4-6), we observe an overall decrease in Re[F(K)] with the increase in frequency, 
i.e., Re[F{K)] < 1 for all k's. 

In Fig. 7 we have also studied the behavior of the dimensionless dynamic permeability K*(lu^) as a function of 
UJ■^,, the dimensionless frequency defined as w* = imW, for the case when a = 0.1. Since a < Uc, we observe sharp 
resonances in the dynamic permability at certain frequences, which can be explained by the non-Newtonian behavior 
of the fluid. 



III. DISCUSSION 

In this paper we have studied non-Newtonian effects on the dynamies of fluids in porous media by calculating a 
flow of Maxwell fluid in a tube, which oscillates longitudinally and is subject to an oscillatory pressure gradient. 
The present study investigates novelties brought about into the classic Blot's theory [1] of propagation of elastic 
waves in a fluid-saturated porous solid by inclusion of non-Newtonian effects. We have used time Fourier transform 
and converted the problem to the frequency domain. We have calculated the dynamic permeability, thus modifying 
the work of del Rio, Lopez de Haro, Whitaker [4] by inclusion of the effect of longitunally oscillating tube walls. 
We investigated how the function F(k), which measures the deviation from Poiseuille flow friction as a function of 
frequency parameter k, is modified by non-Newtonian effects. Our work thus provides a combined theory of non- 
Newtonian flow in a longitudinally oscillating tube, which constitutes the basis for a realistic model of the effects of 
elastic waves in a fluid-saturated porous space. The present analisys clearly demonstrates the existence of a transition 
from a dissipative to elastic regime (as a decreases), in which sharp enhancements of flow (resonances) occur. 

The importance of the current work is two-fold: 

(A) We studied modifications brought about by non-Newtonian effects into Blot's theory. The investigation of the 
function F{k) is important for a number of applications, since F{k) uniquely determines the response of a realistic 
fluid-saturated porous medium to the elastic waves. Thus, determination of F{k) for non-Newtonian (Maxwell) fluid 
is necessary to guide e.g., oil-field exploration applications. 

(B) As we have seen from Fig. 7, non-Newtonian effects cause substantial enhancements in dynamic permeability. 
We were not able to determine in the literature what value of a a natural crude oil would have. However, as it can be 
seen from Fig. 7 for a = 0.1, we can get up to 60 times increase in permeability at certain resonant frequences. Lower 
a's yield even more drastic emhancements. At any rate, we obtained an analytical expression for F{k) (Eq.(9)), which 
can provide the behavior of this function for any given a. 

In the introduction section, the practical impact of the possibility of acoustic stimulation of oil reservoirs has been 
outlined. This result clearly demonstrates that in the crude oil, that can be modelled as a Maxwell fluid, there are 
certain resonant frequencies at which oil production can be increased significantly if the well is irradiated with elastic 
waves at these frequencies. 
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Figure captions 

Fig. 1 Behavior of Re[F(K)] (solid line) and Im[F{K)] (dashed line), as function of k according to Eq.(9). Here, 
a = oo. 

Fig. 2 Same as in Fig. 1 but for a — 10^. 
Fig. 3 Same as in Fig. 1 but for a = 100.0. 
Fig. 4 Same as in Fig. 1 but for a = 10.0. 
Fig. 5 Same as in Fig. 1 but for a ^ 1.0. 



Fig. 6 Same as in Fig. 1 but for a = 10 ^. 

Fig. 7 Behavior of Re[K*{uj^)] = Re[K {oj ;) / K {Q)] (solid line) and Im[K*{uj^)] = Im[K{uj^)/ K{Q)] (dashed line), 
as a function of a;* according to Eq.(6). Here a = 0.1. 
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